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ABSTRACT 

We present the first near-infrared integral field spectroscopy observations of the galaxy MCG- 
6-30-15. The //-band data studied in this paper cover the central 500 pc of the galaxy at the 
best resolution (0".l) so far. The spectra of the innermost regions are dominated by broad 
brackett series emission lines and non-stellar continuum, under which we are able to trace the 
distribution and kinematics of the stars and also the [Fe II] line emission. We find that there 
is a counter-rotating stellar core extending out to 125 pc, which appears to be associated with 
the [Fe II] emission. Based on the mass-to-light ratio, and the presence of this emission line, 
we estimate the age of the central stellar population to be of order of 65 Myr. We show that the 
gas needed to fuel the black hole is, at most, only 1 per cent of that needed to form these stars. 
We derive independent constraints on the black hole mass using the dynamical information 
and determine an upper limit for the black hole mass, Mbh < 6 x 10 7 M q , that is consistent 
with other estimates. 

Key words: galaxies: nuclei - galaxies: active - black hole physics - galaxies: individual: 
MCG-6-30-15 - infrared: galaxies 



1 INTRODUCTION 

The galaxy MCG-6-30-15 is a n elongated lenticul ar (SO) galaxy 
and classified as a Seyfert 1.2 dBoisson et al.ll2002l) . It presents a 
dust lane south of the nucleus and parallel to the photometri c ma- 
. jor axis of the galaxy dMalkan et al.lll998l : Ferruit et al.ll2000h . The 
spectral energy distribution (SED) peaks in the mid-IR and is prob- 
ably due to thermal emission from warm/hot dust grains. In the op- 
tical, the spectrum is dominated by broad Balmer lines and narrow 
forbidden oxygen lines (including [O III] A5007 A ), and the nar- 
row l ine region appears to have a large spatial extent dBennert et al.l 
2006). Stron g internal reddening is observed (E(B — V) ~ 
(0.5 - 1.0) [Reynolds et all 1 19971 : iBoisson etal] |2002h and the 



spect r um show signs of complex obscuration ( Reynolds et al.l 
1 19971 : lBallantvneetai]|2003l) . IRevnolds et ail d 19971) do a multi- 
wavelength analysis of this galaxy, with optical, UV, IR and X- 
ray data and find evidence of a dusty warm absorber (i.e. a col- 
umn of dusty ionised material, usua lly in the form of an out- 
flowing wind) along our line-of-sight. iBallantvne et alj d2003l) ar- 
gue that the dust lane could be responsible for the reddening ob- 
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served and for part of the warm absorber. MCG-6-30-15 has an 
Active Galactic Nucleus (AGN) which has been studied exten- 
sively in the X-ray band (L x (2 - lOkeV ) ~ 4 x 10 42 erg s" 1 
IWinter et al1l2009l: IVasudevan et al.ll2009l) . The non-stellar contin- 
uum can also be observed in the infrared dOlivaetalll 1999b and 
is detected but only marg inally resolved in the radio dNagar et al.l 
ll999l : lM"undell et alj2009h. High angular resoluti on observations in 
the mid-IR drlorst et alj2009l : lGandhi et ail2 009). reaching resolu- 
tions of the order of ~ 0".35, show the nucleus as a point source. 

This galaxy is best-known for having provided the first de- 
tection of a relativistic broadened Fe K a emissi on line gen- 
erated in the inner regions of an accretion disc dTanaka et al.1 
1 19951) . hence providing evidence for the presence of a supermas- 
sive black hole. The X-ray spectral v a riability has been stud- 
ied in detai l (e.g . IRevnolds et ail 1 19951: iNowak & Chian3 200C : 
lUttlev et alj|200a IVaughan & Fabianll2004 iMcHardv et alj|2005 : 



lEmmanoulopoulos et al.ll201ll) and has provided approximate con- 
straints for the mass of the black hole, assuming that the break fre- 
quency scales with the black hole mass. The black hole spin has 
been measured f rom the spectral fitting of the Fe K a emission line 
to be a = 0.989 (Br enneman & Revnoldsll2006l) . i.e. a rapidly ro- 
tating black hole. 

Despite the detailed studies mentioned above, the mass of the 
black hole responsible for the AGN activity is still not very well 
constrained. The galaxy is too distant for current instruments to 
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Figure 1. HST and SINFONI images of MCG-6-30-15. Left: WFC2 broadband image of MCG-6-30-15 from lMalkan et alJI 19981 . Right: SINFONI //-band 
data cube collapsed to two dimensions for plotting purposes. The SINFONI field-of-view is rotated by 25 degrees clockwise relative to the HST image. 



resolve the black hole sphere of influence. Recent work to mea- 
sure the black hole mass of MCG-6-30-15 has been done by 
McH ardvetal.l l l2005h who discuss several methods (Mbh — o re- 
lation, virial type relations and spectral variability) to determine 
the mass. They measure the stellar velocity dispersion from the Ca 
II absorption lines up to a radius of R e /8 and obtain a value of 
cr = 93.5 km s . Us i ng the various M bh — cr relations available 
dTremaine et~alll2002l : If errares el Eooi and constraints from vari- 
ability, the mass is fo und to be in the rang e Mbh ~ (3 — 6) x 10 6 
M0. Another study feennert et al.ll200q) . finds somewhat larger 
black hole masses (0.8 - 3) x 10 7 M Q , after using several meth- 
ods to determi ne it (virial determinations and Mbh-c relation). 
More recently, Iwinter et al.l d2009b used the met h od ba sed on the 
Mbh— Lk relation outlined in Mushotzkv et al.l d2008l) to calcu- 
late a black hole mass of Mbh = 2.3 x 10 7 MM. fVasudevan et al.l 
(2009) refined this method and determined a similar black hole 
mass: Mbh = 1.8 x 1O 7 M . 

In terms of the stellar population, iBoisson et al.l d2004l) finds 
the stellar content to be d ominated by an old population. From 
analysis in the ultra-violet, Bon atto et alj (J2000) find a dominant 
old bulge stellar population (~ 10 Gyr) with indications of a se- 
ries of previous bursts of star formation distributed in age among 
the younger stellar p opulation (< 500 Myr). Optical analysis by 
iBennertet all [2006 ) reach the same conclusion. 

The advent of integral field spectroscopy (IFS) has opened 
a new perspective on the black hole fuelling process and 
on the relation between the stellar properties of the galaxies 
and the AGN activity. Studies using this technique have un- 
veiled features in the gas kinematics which are likely due to 
the dynamics of inflowing and outflowing gas. Observations 
of local Seyfert galaxies show evidence for gas inflowing or 
outflowing on a scale of tens to hundreds of parsecs (e.g.: 
Prieto et al.ll2005L iFathi et al.ll2006l IStorchi- Bergman n et alj|2007l 



Mtille r~anchez et alJl201l(K Combined with adaptive optics, IFS 
provides the opportunity to map with high resolution the in- 
ner regions of nearby galaxies and determine the ma ss of the 
central black hole (e.g. NGC3227 - iDavies et al.1 I20061. Centau- 
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iNowak et al.l2008l NGC524 andNGC2549 - iKrainovic et alj2009l 



M87 - iGebhardt et alj201 lh. We can also emp loy IFS to investigate 
the AGN nature (e.g. lValencia-S. et al.|[2012l) and th e relationship 
betwe en the star formation and the AGN activity (e.g. lDavies et al.l 
120071) . The high level of detail and amount of information that 
comes out of these observations makes them a very useful probe 
of the physics in the central regions of galaxies. 

MCG-6-30-15 is an ideal system in which to investigate the 
physics of black hole accretion, but there are still many aspects 
about its nature that remain unknown. By observing the inner re- 
gions of this galaxy, we can study its stellar and dynamical proper- 
ties and draw conclusions on the black hole and AGN activity. In 
this work we present the first study of this galaxy using IFS data in 
the //-band, and determine the kinematic properties in the central 
regions of MCG-6-30-15. Section [2] gives the details on the data 
reduction. In Section|3]we present the results obtained for the stel- 
lar kinematics and gas dynamics traced by the [Fe II] emission line. 
We use dynamical arguments to constrain the black hole mass and 
the star formation history of the galaxy. The results are summarised 
in Section[4] 

We adopt the standard cosmological parameters of Ho = 70 
km s" 1 Mpc -1 , fi m = 0.27 and D. A = 0.73. The scale for 
our measure ments is approxim ately 0.158 kpc/" at a redshift 
2 = 0.0077 dFisheretal.lll995l value quoted at the NASA/IPAC 
Extragalactic database ). The distanc e to MCG-6-30-15 used in 
this paper is 33.2 Mpc jWrightl2006T) . 



2 DATA ANALYSIS 
2.1 Data reduction 

We obtained spectroscopy data for our source, wi th the near-IR in- 
tegral field spectrograph SINFONI on the VLT (Eisenhau er et al.l 
2003i; iBonnet et al.l f2004). using a natural guide star for the adap- 
tive optics (AO). The observations were carried out in the //-band 
mainly to avoid the contamination from the non-stellar contin- 
uum. In the near-infrared, the relative contribution of the non-stellar 
emission compared with the stellar emission is lower than in other 
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Figure 2. Spectral data integrated in the central 6x4 pixel region of the 
field-of-view and the best fit spectral components as a function of the rest- 
frame wavelength. Top: Data (black) and best fit model (red). Bottom: Spec- 
tral components used in the fit, broad Brackett emission lines (black) and 
fitted sum of stellar templates tracing the absorption features (green). The 
scale is normalised to the continuum value. The excess at around 1.55 fim 
corresponds to the position of another broad line in the Brackett series. We 
tried to extend the fit to this line with a sixth Brackett emission line but due 
to its low intensity the fit was not very good. We decided then not to fit this 
line and exclude its wavelength range from the subsequent data analysis. 



wavelengths and less affected by extinction. The //-band, in partic- 
ular, sho ws less nuclear emi ssion reprocessed by dust than in the 
iT-band dBoisson et alJl2002h . The //-band spectra with R~3000, 
were taken on three nights, (April 2006, 4th, 5th and 20th), with a 
pixel scale 0".05"x0".l, sampling a field-of-view of 3"x3" (right 
panel of Fig. [T). Each exposure was of 300 seconds, alternating be- 
tween the object (A) and the sky (B) in the sequence ABBA. This 
sequence was repeated eight times. The total on-source exposure 
time is lh20m. 

The data were reduced using the ESO pipeline for SIN- 
FONI version 3.8.3, and independent routines when necessary. 
Using the pipeline and the calibration files, we obtained the bad 
pixel maps (non-linear pixels, hot pixels and reference bad pixels), 
master dark and master flat frames (recipes SINFO_REC_DETLIN, 
SINFO_REC_MDARK and SINFO_REC_MFLAT). These routines fail 
to identify some of the bad pixels, so we visually inspected our raw 
frames and manually selected extra bad pixels, in particular hori- 
zontal bad pixel lines that can be seen in the first pointings of each 
of the four ABBAABBA sets of observations. 

We subtract the master dark from the raw object and sky 
frames, and use the recipe SINFO_REC_JITTER to correct for 
bad-pixels and distortions, flat-field the data, wavelength cali- 
brate and reconstruct 32 data cubes for object and sky. The 
SINFO_REC JITTER routine is run with the sky subtraction options 
turned off, to allow us to correct the obje ct data cubes for sky 
emission following the method described in IPaviesI (120071) . From 
this routine we obtain each input data cube corrected for the sky 
emission. Our final 16 object data cubes are then checked for any 
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Figure 3. Equivalent width of the CO (6-3) 1.6187 fim absorption line 
measured in integrated annuli at distance r from the AGN position. The 
change in equivalent width is due to non-stellar continuum associated with 
the AGN. 



remai ning bad pixels usi ng a 3D Laplacian edge detection IDL 
code ( Davies et ajJ 120101) . which is a 3D version of LA cosmic 
( Ivan DokkurrfeOOlh . 

In the near-IR, the atmosphere greatly affects the transmission 
and there are several telluric absorption signatures visible in the 
spectrum. We observed four telluric standard stars in the //-band 
but we found those were not suitable due to low signal-to-noise 
ratios. Instead, we divide our data cubes by a normalised theo- 
retical template telluric spectrum re-scaled to the air-mass of our 
observations. The last step was to combine the 16 individual ex- 
posure cubes into a final cube. We re-sample each cube to a half- 
pixel scale both in x and y dimensions and offset them using the 
values recorded in their headers (CUMOFFSET and TEL TARG AL- 
PHA/DELTA), verifying that the bright galactic nucleus is aligned in 
every spatially shifted cube, and then median combine them. The 
final data cube (collapsed to two dimensions for plotting purposes) 
is shown in the right panel of Fig. Q] The pixel scale of the final 
reduced cube is 0".05 x 0".05 and the horizontal x-axis is aligned 
with the galaxy's major axis, which means a tilt of ~ 25° in rela- 
tion to the East- West direction. 



2.2 Flux calibration 

The flux calibration of our data cube is done based on the 2MASS 
//-band image of the galaxy, since our standard stars were not suit- 
able (low S/N or saturated). We integrate the counts in a 3" diame- 
ter circular aperture in the 2MASS image, which gives an observed 
magnitude of mg = 11.3 and a flux of 3.37 x 10~ 14 W m~ 2 /zm -1 
(Cohen et al]|2003h . We then integrate the spectra in our data cube 
in a 3" aperture centred at the same physical position as for the 
2MASS image. We compare the median count number of our in- 
tegrated spectrum in the range A = [1.509 — 1.799]/im with the 
flux from the 2MASS image and use this as our count to flux cal- 
ibration. Our calibration was checked using different apertures in 
2MASS and is consistent to ~ 20 per cent. 



2.3 Emission components 

In the case of MCG-6-30-15, the broad hydrogen Brackett AGN 
emission lines, vfwhm ~ 1800 km/s, are very strong in the 
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Figure 4. Panel with decomposed emission. Left: Stellar continuum. Centre: Broad Hydrogen Brackett emission lines. Right: Non-stellar continuum. Stellar 
continuum contours are overlapped on the three panels and provide a comparison between the stellar distribution and the unresolved AGN components (broad 
Brackett emission lines and non-stellar continuum). 



central pixels and dominate over the wavelength range of inter- 
est, making the stellar absoiption features hard to identify. The 
following method is adopted to remove the broad emission lines. 
We first select an integrated central 6 by 4 pixel region where the 
broad lines are stronger and fit the five strongest emission lines 
with Lorentzian profiles, since these provide a good fit to the lines 
dVeron-Cettv et all l200lh . We fix the atomic parameters and the 
velocity width of the line (relative half-width at half-maximum and 
relati ve intensity - using as an approximation Hum mer & Storevl 
1987 for case B recombination). To improve the fitting we also in- 
clude a second order polynomial function to fit the non-stellar con- 
tinuum and a giant-star template spectrum convolved with a Gaus- 
sian as a basic model to the stellar continuum and its absorption 
features. The stellar template is of a Ml III star observed previ- 
ously with SINFONI in the 100 mas scale. It was chosen due to 
its deep CO absorption lines which are also observed in our spec- 
tra. The plot in Fig. [2] shows the spectral components fitted to the 
data. The best-fit parameters for the Brackett lines obtained from 
the integrated spectrum are used as initial guesses when fitting ev- 
ery spaxel. The spaxel-by-spaxel fitting determines the line param- 
eters for the Brackett emission at each spatial position and removes 
it from the data cube. 

With this approach we can decompose the emission into the 
individual AGN and stellar contributions. The broad Brackett emis- 
sion is due to the AGN, while the absorption lines are due to the 
stellar population. However, the continuum is composed of stellar 
emission and AGN heated dust. To separate these two contributions 
we use the equivalent width of one of the strongest absorption lines 
present in our data: CO (6-3) A = 1.6187. The argument is that the 
deviations of the equivalent width from the intrinsic value will be 
caused by dilution due to the non-stellar continuum. By measuring 
the equivalent width at different radial positions and comparing the 
measurements to an intrinsic value, it is possible to determine the 
relative AGN and stellar contributions at each position. The stellar 
fraction will be, 

/stellar W = ^^-(r), (1) 
vVintr 

where Wohs is the observed equivalent width and Wintr is the in- 
trinsic equivalent width. We measure the equivalent width of the 
line in integrated radial bins, and plot it as a function of distance 
from the nucleus in Fig. [3] We fit a Gaussian to the distribution and 



use it to decompose the continuum. The maximum equivalent width 
(~ 3.1A at r= 0".8) is taken as the intrinsic value for our case. 
Although it is lowe r than what has been ob served before (~ 4A) 
iDavies et ai] fc007h . lValencia-S. et all ( |2012|) . this could be due to 
the presence of dust fa rther out in the ga laxy. Lower values have 
also been observed by lOliva et ail d 19991) using slit-spectroscopy 
on MCG-6-30-15. In Fig. [4] we show the decomposition maps of 
the observed emission into the separate contributions from the stel- 
lar continuum (left), the Brackett broad emission (centre) and the 
non-stellar continuum (right). In Fig.|5]we show the cuts along the 
major axis, passing through the AGN position determined by the 
peak of the Brackett emission (offset = arcsec). From the plot, it 
is clear that the Brackett emission and the AGN continuum have 
approximately the same distribution, which is what we would ex- 
pect since they are both spatially unresolved l lBennert et al 1 120061 
estimate the broad line region to be 9 light-days across). The stel- 
lar continuum shows more extended wings which correspond to the 
stellar population spatially distributed in the galaxy. 

2.4 PSF and spectral resolution 

We determine the point spread function (PSF) from the Brackett 
emission map. The broad Brackett emission is coming from a very 
small region around the AGN and it is not resolved in our data. 
The spatial extent of the observed emission will be caused by the 
instrument and the atmosphere and we can measure the PSF by 
modelling its distribution. In our data the PSF is well modelled by 
a double Gaussian: a Gaussian with small width to model the peak 
(FWHM = 1 .4 pixels ~ 0" .07) and a wider one with (FWHM = 4.0 
pixels ~ 0".2) to model the wings. The non-stellar continuum can 
also be used to determine the PSF since the dust-emitting region is 
unresolved. We obtain similar results: FWHM = 1.4pixels~ 0".07 
and FWHM = 3.7 pixels ~ 0".18. 

The spectral resolution is obtained from the sky lines. These 
emission lines are present in our sky exposures, and their width 
is due to SINFONI's instrumental broadening. We use the com- 
bined sky cube and fit an unblended sky line with a Gaussian. The 
half width at half maximum (HWHM) is taken as the instrumental 
broadening (ajnst = HWHM / ^/21n(2)). We repeat this procedure 
for all the spatial positions in a masked cube, to exclude the noisy 
outer regions of the field-of-view. The mean value for the instru- 
mental broadening is <Ti na t = 60 km/s (3.2 A at A ~ 1.6187/im), 
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Figure 5. Distribution along the major axis of the galaxy. The zero offset 
corresponds to the AGN position. Top: Brackett emission. Middle: Non- 
stellar continuum. Bottom: Stellar continuum. The stellar contribution has 
a broader core and clearly extends to a larger radii than the one from the 
AGN. The orientation is as for Fig. [4] 



with variations between 55 km/s and 67 km/s. The unblended line 
is selected in a wavelength region around the central wavelength of 
[Fe II]. However, in the wavelength range of our analysis, we note 
that the maximum instrumental broadening variation is 5 km/s. 



library contains only ten G, K and M giant stars, but as we found 
previously, K and M stars reproduce our spectra well and are suit- 
able to model the stellar absoiption features. To remove the effects 
of differential instrumental broadening between the SINFONI data 
and the templates, we broaden the stellar templates (FWHM tcmp 
= 8 pixels ~ 3.2 x 10~ 4 /im) to the resolution of SINFONI by 
convolving them with a Gaussian kernel (PSF_GAUSSIAN in IDL) 

of cr conV oi = y'ofnst ~~ °"tem p - °fnst is tne SINFONI instrumental 
resolution which we measure from the sky emission lines. 



2.5.1 Binning and masking 

We use the strength of the absoiption lines to bin the data spatially, 
and in this way avoid being dominated by the noise from the AGN 
continuum. Usually the signal-to-noise is quoted as the 1<t noise 
value in relation to the stellar continuum or total continuum. Here 
we define our signal-to-noise ratio as the depth of the absorption 
lines divided by the la noise level, since this is ultimately the mea- 
surement we are i nterested in. We bin our data using the 2D Voronoi 
binning routine o f lCappellari & Copinl J2003h to an average signal- 
to-noise S/N = 5 (between the amplitude of the CO (6-3) 1.6187 
absorption line and the noise RMS in a line free region of the 
continuum). This corresponds to a S/N ~ 50 when measured in 
relation to the stellar continuum level. A more accurate signal-to- 
noise measurement will be determined in the next sections from the 
residuals of our fit. 

Individual pixels in the regions where the flux decreases to 
1/125 of its peak value were masked out and excluded from the 
analysis. These pixels in general correspond to the outer regions 
of our field-of-view. We also masked out a ~ 6 pixel wide region 
parallel to the major axis at around 0".8 from the centre of the field- 
of-view, which corresponds to an illumination artefact in SINFONI. 



2.5 Extracting stellar kinematics 

The velocity map and stellar velocity dispersion distribution are 
determined from our fully reduced and br oad line subtracted data 
cube, by fitting the spectra with p PXF (jCappellari & E mselleniil 
120041 : Ivan der Marel & Franxlll993l) . This IDL routine uses a set 
of stellar templates to find the best weighted spectral combination 
and line-of-sight velocity distribution to fit the input spectrum. The 
line-of-sight velocity distribution is based on the form of a Gauss- 
Hermite series expansion, with the first two moments being the ve- 
locity offset and the velocity dispersion. It is also possible to fit 
higher moments such as the asymmetric and symmetric deviations 
from a Gaussian (namely the h3 and /i4 parameters which are re- 
lated to the skewness and kurtosis of the distribution). In this work 
we fit only the first two moments of the distribution due to the lim- 
ited signal-to-noise of our data. 

In the //-band, there are not many medium or high resolution 
stellar templates. Our data have a resolution R ~ 3000, which is 
higher than most of the templates found in the literature. For this 
reason we use different sets of stars to determine the stellar pop- 
ulation and t he lin e-of-sight distribution. The stellar templates of 
iMever et alj l ll998l) . which include several spectral types and lumi- 
nosity classes, have spectral resolution of the order of our data res- 
olution (R ~ 3000), but we encountered problems when determin- 
ing the velocity properties. We use this library to determine that our 
spectra can be well fit based on a set of K and M stars. To determine 
the velocity properties we use the more rece nt stellar templ ates at 
higher resolution of R ~ 5000 obtained bv lLe et alj feoill) . This 



3 DISCUSSION 

3.1 Stellar kinematics 

The stellar kinematic properties are determined by running pPXF 
on our binned da ta cube. The input stellar library is the one men- 
tioned above ihe et alj|201lh . and we include an additive poly- 
nomial of fourth order to model the AGN continuum. The emis- 
sion lines of [Fe II] and a residual from the telluric subtraction are 
masked out when fitting the spectra. The code output includes the 
best fit spectrum at each pixel, the velocity offset and the velocity 
dispersion. We do not fit the higher moments h3 and h4 as men- 
tioned in the previous section due to our low signal-to-noise. 

We first run pPXF in the integrated field-of-view, excluding 
the central regions dominated by the AGN continuum to increase 
the signal-to-noise ratio. We plot the results, which correspond to 
our maximum signal-to-noise ratio for these data (S/N ~ 10 for the 
deepest absoiption features or S/N ~ 150 in relation to the total 
continuum) (Fig. [6}. Our stellar templates provide a good fit to the 
data, with the main contributions of the stellar templates of a K5 III 
and a M0 III star. We removed an absorption feature that the stellar 
templates were not able to fit, it coincides with the wavelength of 
the Fe I stellar absorption and a telluric feature. The fact that we are 
not able to fit it could be due to telluric residuals and/or a metallicity 
higher than the metallicity of our templates (which do not include 
super-metallicity stars). 

In Fig. [7] we show the spectral data and spectral fitting at 
four different regions of the field-of-view (numbered boxes in the 
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Figure 6. Spectrum and fitting results versus rest-frame wavelength in (im. In black, integrated spectrum over a radius of ~ l".l excluding the central 6x4 
pixels. In red, pPXF best fit combination of stellar templates and 4th order polynomial (to model the AGN continuum). The residuals are shown in grey and 
have been offset from zero for plotting purposes. The vertical shaded regions indicate wavebands which were excluded from the fit, including the [Fe II] 
emission lines at A = 1.644^tm and A = 1.677/xm and the absorption feature at A ~ 1.654/^m (see text for details). 



top panel). The spectra plotted in the bottom panels, show the in- 
tegrated spectrum (black), best-fit model (red) and the residuals 
(grey), for each of the four different regions. It is clear that the stel- 
lar features are deeper close to the major axis of the galaxy (regions 
1 and 3), where we also see the stronger stellar continuum. 

We follow the same procedure as described above but now 
across the entire field-of-view on a spaxel-by-spaxel basis, fitting 
the wavelength range A rcst (1.57 - 1.716) ^m, to obtain a map 
of stellar velocity and velocity dispersion. We show the results in 
Fig. [8] The systemic velo city is determined usin g the method de- 
scribed in Appendix C of Krain ovTc et al.l d2006l) . using the veloc- 
ity maps of the central r > 1". We obtain best fit values of V sys = 
2396 km/s and a kinematic PA of 1 12± 16 degrees. The velocity is 
corrected for the g alaxy's sy s temi c velocity, but not for the inclina- 
tion. F rom work of Lauberts ( 1982) compiled in lde Zotti & GaskeiH 
( 1985), the axis ratio for this galaxy is b/a — 0.58. The axis ratio 
of our stellar continuum is consistent with the value quoted above, 
indicating that the axis ratio at small and large scales is the same. 
From the axis ratio of the stellar continuum (and assuming a flat 
geometry), we derive an inclination of ~ 55 degrees. Our rotation 
velocities could be up to 49 km/s, instead of the observed ~ 40 
km/s. As described in Section |231 the velocity dispersion obtained 
already takes into account the instrumental broadening. In the cen- 
tral regions (r < 0".l) around the black hole, the non-stellar con- 
tinuum is very strong, and the relative intensity of the stellar ab- 
sorption lines is low, which increases the error in determining the 
velocity dispersion. For this reason we masked out the inner r < 
0".l (white bins in Fig. [8}. The errors in the parameters are mea- 
sured using a Monte Carlo approach. We generate a random set of 
100 spectra with the same noise properties (S/N ~ 5) as our galaxy 
spectra and fit them with pPXF. The starting values and the wave- 
length limits of the fit are changed randomly as well. The final la 
deviations in each parameter are taken as our errors for a typical 
spectrum with S/N ~ 5 (which corresponds to the S/N in each bin). 
We find errors of ±4.5 km/s for the mean line-of-sight velocity and 
±5.1 km/s for the velocity dispersion. 



The counter-rotating core 

The results show a low rotational velocity when compared with the 
velocity dispersion. The zero point in the velocity map corresponds 
to the systemic velocity of the galaxy, negative velocities are as- 
sociated with blue-shifted absorption lines and positive velocities 
with redshifted absorption lines. The central ~ 1.4 arcsec radius 
region shows the evidence of two distinct kinematic components. 
The inner and outer regions of our map indicate two different and 
counter-rotating components, with the r < 0".8 ~ 125 pc region 
rotating clockwise and the r > 0".8 rotating counter-clockwise. 
The two kinematic components have similar absolute velocity val- 
ues (40 — 50 km/s). The rotation of the counter rotating core is 
seen clearly in the figure. The outer region is only beginning at the 
limits of our analysed region of the field-of-view - corresponding 
to the region above the flux threshold. It is visible in the top right 
(yellow/red zone) and top left (blue zone) of the velocity map. The 
velocity dispersion map presents values that decrease inwards and 
are flat at the centre. It has values of ~ 130 km/s in the outer re- 
gions r > 0".8 but drops by ~ 30 km/s in the inner regions, which 
coincides with the counter-rotating core. The dispersion values are 
in general higher than the line-of-sight velocity, which suggests that 
the galaxy core is not rotationally supported. 

Although we do not use the h3 distribution in our analysis 
due to the low signal-to-noise ratio, we did a test to check if it was 
consistent with the velocity map of the top panel of Fig. [8] A value 
of h3 different from zero indicates that the line-of-sight velocity 
distribution profile deviates from a Gaussian distribution. The value 
of h3 is related with the skewness, and measures the asymmetry 
in the distribution. As expected, the h3 map is very clearly anti- 
correlated with the line-of-sight velocity distribution for r < 0".8 
and in the upper regions of the map where the velocity changes 
sign (although not as clearly due to the higher noise level in those 
pixels). 

The general kinematic properties mentioned above, indicate 
that we are most likely in the presence of a kinematically distinct 
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Figure 7. Integrated spectra and best fit model for four different regions of the field-of-view. The size of each region is 5 x5 pixels indicated with the green 
boxes in the top image. The bottom four panels show the spectra (black), the best fit model using pPXF (red) and the residuals (grey) offset for plotting 
purposes. The vertical shaded bars indicate regions excluded from the fit. 



core. These type of systems are fairly common in ellipticals but 
cou nter-rotating systems are rarer ( < 10 per cen t ) in SO type galax - 
ies dKuiiken et ail 1 199d : iKrainovic et allfcoi ll ; iBois et alj 1201 ll) . 
The drop in the velocity dispersion occurs for r < 0".8 and is 
possibly related with the counter-rotating core. This phenomenon 
has been observed for another SO galaxy with a kinematically- 
distinct core, NGC 7332, where the velocity dispersion increases 
from r ~ 10" or 1 kpc towards the centre of the galaxy, but 
drops by ~ 10 km/s when in th e region of the counter-rotating 
core dFalcon-Barroso etal 1 120041) . Our data cover a smaller radii 
(r max ~ 1".4 or 200 pc) than the work of iFalcon-Barroso et al.l 
(2004), which allow us to measure the drop in velocity dispersion 
in the inner regions, but does not permit us to evaluate if the ve- 
locity dispersion decreases again as we move to larger radii. The 
dynamically decoupled core could be the result of an inflow of gas 
into the central regions of the galaxy via, for example, a past minor 



merger event. For our target, the majority of the stellar population 
is old and fairly homogene ous, not showing any variations with 
radius I TBoisson et "ai] |2002l) . There is nevertheless indication of a 
series of previous star formation bursts dBonatto et al.ll200"oh . The 
evolution for this galaxy is not expected to be through major merg- 
ers, since only about 2 — 2.5 per cent o f SO galaxies are e xpected 
to have had a major merger in their past (IParrv et al.l2009l) . On the 
other hand, secular evolution is expected to dominate the evolu- 
tion of Narrow Line Seyfer t 1 galaxies such as MCG-6-30-15 (e.g. 
lOrban de Xivrv et al.ll201 lh . 



3.2 Gas kinematics 

When subtracting the stellar and AGN continuum from our data 
cube using the results from pPXF, we are left with the gas emis- 
sion spectra at each spatial location. The A = 1.644 /im forbidden 
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Figure 8. Two-Dimensional map of line of sight velocity and velocity dispersion from the spectral fitting with pPXF. The orientation is the same as in the 
right panel of Fig. 1. The white bins were masked out due to the high non-stellar contribution and the black cross indicates the AGN position. Top: Line of 
sight velocity map. Positive and negative velocities correspond to movement away and towards the observer respectively. The values were corrected for the 
systemic velocity of the galaxy but not for the inclination (~ 55 degrees). To note that the values are more uncertain in the central pixels due to the presence 
of a high AGN continuum. In the central regions the stars rotate in the clockwise direction while in the outer regions they rotate in the opposite direction (seen 
as a change of velocity direction in the top left and top right regions of the map). Bottom: Map of velocity dispersion. The velocity dispersion is higher in the 
outer regions, decreasing when approaching the nuclear region. The higher values of dispersion occur in a plane, showing little vertical extension. The values 
of velocity dispersion have been corrected for the instrumental broadening. 



emission line of [Fe II] is the strongest emission line, and clearly 
observed by eye in the spectra. In Fig.[9]we plot the fit to the emis- 
sion line integrated in two regions of the same size (10 x 10 pixels) 
offset by 0".5 from the AGN position on the East and on the West 
side of the nucleus. As we can see from the plot, the intensity is 
larger on the West than on the East side. The line velocity is also 
different, the West side is blueshifted (—8 km/s) and the East side 
is redshifted (+129 km/s). The lines show a broadening of 49 km/s 
and 74 km/s on the East and West side of the nucleus respectively. 
In the individual pixels, the S/N of the line is not as strong. To get 
the spatial distribution of [Fe II], we fit a Gaussian to this line and 
fit the continuum usin g regions on the left and righ t of the line using 
the algorithm MPFIT (Mar kwarddl2009l,lMorelll978l) . We require a 
signal-to-noise ratio S/N > 3 in relation to the RMS scatter of the 
spectrum, a condition that in general holds in the central regions 
of the galaxy and coincides with the kinematically-distinct core. 
We show the two-dimensional line properties in the three panels of 
Fig. HO] the top panel shows the velocity offset, the middle panel 



the velocity dispersion and the bottom panel the line intensity. In 
these plots, the velocities are measured in relation to the systemic 
velocity of the galaxy. The instrumental broadening is subtracted in 
quadrature using the dispersion of an unblended sky emission line 
close to the [Fe II] line emission. We excluded regions that showed 
velocity dispersion values similar to the instrumental broadening. 



The lines show a gradient from the minimum values 
blueshifted by ~ — 30km/s on the West side of the nucleus to 
the larger redshifted values (+120 km/s) on the East side. We ob- 
serve an asymmetry in the intensity map, with larger values on the 
West side of the nucleus than on the East. The velocity dispersion 
is harder to interpret, it is higher on the West side of the nucleus 
(~ 60 km/s) but decreases around the AGN position. It seems to 
increase again on the East side although the signal-to-noise ratio 
there is not as good. 
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Figure 9. [Fe II] emission integrated in two regions (10x10 pixel each) 
of the field-of-view and plotted as a function of the rest-frame wavelength. 
East (top panel) and West (bottom panel) of the nucleus, along the galaxy 
major axis. The gas shows a different velocity in both regions, as can be 
seen from the wavelength shift in the plot. The integrated Western region 
shows a larger velocity dispersion and intensity than the emission from the 
Eastern region. The specific values are given in the text. 



The origin of the [Fe II] emission 

Since the signal in the data is weak, we cannot say with certainty 
what causes this [Fe II] emission, but we can discuss the two pos- 
sibilities: an AGN outflow or supernova shocks. The intensity of 
[Fe II] is asymmetric, with stronger emission on the West side of 
the nucleus. The velocity map for [Fe II] shows that the gas has a 
distribution along the main axis of the galaxy and the values show 
a deviation from the rotation velocity of the stars. T he distribution 
of [O III] 5007 A emission us ing HST observations dSchmitt et al.l 
120031 and iBennert et al . 2006) show an elongation along the ma- 
jor axis of the galaxy, as we see in the [Fe II] distribution. The [O 
III] emission also shows an asymmetric distribution but on a larger 
scale than the field-of-view of our observations. The fact that we see 
a velocity structure distinct from the stellar rotation pattern, and an 
elongation similar to the [O III] emission, could indicate that we 
are in the presence of an AGN driven outflo w, with [Fe II] and [O 
III] em ission tracing the same geometry (e.g. Mullcr -Sanch ez et al.l 
l]). The gas kinematics could be consistent with emission from 
a cone that is outflowing away from us. 

There are nevertheless some caveats in this hypothesis. The 
velocity offset from the rotational velocity is not very large, ~ 10 
km/s on the West side and ~ 80 km/s on the East side of the nu- 
cleus (where the S/N is lower), and could be due to uncertainties in 
the velocity measurements. The transition in velocity as it moves 
across the nucleus is fairly smooth and resembles a rotation pat- 
tern. Our data samples a much smaller spatial scale than the [O 
III] maps, and we cannot check if the asymmetry of the [Fe II] 
intensity is observed at the larger scales as well. We are limited 
to the small scale of our observations and the low signal-to-noise. 
With new data on a larger scale, it would be possible to investi- 
gate if the asymmetry in the gas emission holds for larger radii. 
It is also known that at the very small scales probed by the X- 
ray observations of this source, there are spectroscopic signs of a 
warm absorber (e.g. iRevnolds et al]|l997L iBallantvne etafll2003l) . 
and hence a form of outflow. It would be interesting to investigate 
this topic further to determine if there is any connection between 






Figure 10. Panels showing the [Fe II] A = 1.644^tm emission properties 
of MCG-6-30-15 obtained by fitting the line with MPFIT. The white cross 
marks the position of the AGN and the map orientation is the same as for the 
right panel of Fig.[T] The regions referred to in the text as to the East and to 
the West of the nucleus are measured along the major axis. Bins with S/N 
< 3 where excluded and bins for which the velocity dispersion is lower 
or similar to the instrumental broadening were also excluded (the region 
surrounding the black hole position did not obey this condition and hence 
was masked out). Top panel: Velocity corrected for the systemic velocity 
of the galaxy in km/s. Middle panel: Velocity dispersion corrected for the 
instrumental broadening in km/s. Bottom panel: Line intensity. The line 
intensity is higher on the West side of the nucleus than on the East side. The 
velocity shows an asymmetry, with higher absolute values on the East side 
of the nucleus than on the West. The two hypothesis for this [Fe II] emission 
(AGN outflow and most likely supernova induced shocks) are discussed in 
the text. 
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the accretion disc scale outflow and the tens of parsecs of the [Fe 
II] emission. 

The [Fe II] emission is excited by electron collisions and is a 
good tracer of shocks. In the case of MCG-6-30-15, the extended 
emission could be caused by nuclear mass outflo w shocks with am - 
bient clouds, or by supemova-driven shocks dMouri et al.l 2000). 
The spatial distribution of the [Fe II] coincides with the region 
where we see the strongest signal-to-noise in the stellar absorp- 
tion features. This and the arguments presented above, imply that 
most likely the [Fe II] emission is due to shock fronts from su- 
pernova remnants. These shocks destroy the dust grains which al- 
lows for a higher abundance of iron in the gas phase. The iron is 
then ionised by the interstellar medium. The [Fe II] line is there- 
fore a good indicator for the pre sence of supernova remnants and 
a tracer of the supernova ra te dMoorwood & Olival Il988b IColinal 
ll993l : lRosenberg et alj|2012h . 

In this case, the [Fe II] emission can be used to constrain the 
supernova rate and investigate the properties of the stellar popula- 
tion of MCG-6-30-15, as we will discuss in Section[331 



3.3 Central stellar population 

From the observations and results of the previous sections, we can 
investigate the properties of the stellar population of MCG-6-30- 
15. Assuming that the supernova shocks were the main excitation 
mechanism for [Fe II], we can determine the supernova rate (SNR) 
based on the [Fe II] 1.644 /im line flux. We measure a flux (inte- 
grated in the r < 0".8 region) of F[ Fe //] 16 44 



8.2 x 10~ 16 erg 



s" 1 cm -2 which corresponds to a luminosity of I/[Fe//]i.644 = 
1.2 x 10 38 erg s . This [Fe II] is not particularly strong, its flux 
is comparable to the low er values found for narrow lin e Seyfert 1 
galaxies in the sample oflRodrfguez-Ardila~e t al. (2004). Using the 
equations from Rosen berg et al.ld2012 ): 



l og ®™ = ( .89 ± 0.2) x log _ (36.19 ± 0.9) (2) 



[yr~ 



[ergs - 



and t he theoretical ratio [Fe III 1-6 4 /xm/[Fe II] 1.26 pm = 0.7646 
from iNussbaumer & Storevl dl988l) . we calculate a supernova rate 
of SNR = 6.6 x 10 -3 yr" 1 . 



We used the evolutionary synthesis code S TARS dSternbergl 
1 1998MThornlev et al.l200ol:lsternberg et al.l2003l) to follow the evo ; 



lution of a s tellar pop ulation as in 
iDavies et al] 120051 and iDavies et al 



Forster Schreiber et al. I [20031; 



20071 . In the code, the star 



formation declines exponentially with a cha racteristic timescale , 
which we assume to be 10 Myr, adopted from lDavies et al.l d2007t) . 
Assuming a Salpeter Initial Mass Function and this star formation 
timescale, STARS calculates the distribution of stars in function 
of age. The code results are scaled to the case in study by com- 
paring the A'-band luminosity in the model with our observed H- 
band luminosity. The conversion between H and A'-band luminosi- 
ties is done using the magnitude relation H - K = 0.15 which is 
very weakly dependent on the stellar age. The code calculates pa- 
rameters of interest for our investigation, such as the supernova 
rate, A'-band luminosity, mass loss rate and total stellar mass at 
any time. Using the relation between predicted SNR and the stel- 
lar age from the code, we determine the age of the starburst that 
produced the SNR observed. From our SNR we obtain an age of 
6.5 x 10 7 yr. We note that although two values for the age are con- 
sistent with the supernova rate measured, the alternative younger 
age (3 x 10 6 yr) is inconsistent with the presence of deep absorp- 
tion features that indicate the presence of late-type stars. With a 
stellar age of 6.5 x 10 7 yr we obtain from the model a mass-to- 
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Figure 11. Multi-Gaussian Expansion (MGE) modelling of the galaxy sur- 
face brightness. Data (black contours) compared with the MGE model (red 
contours). 



light ratio of M/Lh = 0.98 Mq/Lq. From the model we also 
know the mass loss rate, which can be compared with the mass 
accreted by the black hole. The stronger stellar and [Fe II] emis- 
sion coincide with the counter-rotating core. This type of structure 
is thought to be formed by an inflow of new gas via, for example, 
a minor merger that generates a stellar population with different 
properties in the nucleus of the galaxy. New gas may have gener- 
ated the stellar population we observe and model with STARS, but 
can also have fuelled th e AGN. From th e bolometric luminosity 
(~ 1.3 x 10 44 erg s" 1 ; IVasudevan et alj |2009). and assuming an 
efficiency of e = 0.1 we can obtain an estimate for the mass ac- 
cretion rate M — Lboi/ec 2 . For this AGN the mass accretion rate 
is around 0.022 M Q yr _1 , which gives, assuming a constant mass 
accretion rate during the age of the starburst, a total mass accreted 
of: 



Aft,* = 6.5 x 10 7 x 0.022 = 1.4 x 10 6 M Q 



(3) 



The model-derived total mass used to form stars is 2 x 10 Mq. 
This means that of the new inflow of gas, only around 1 per cent of 
the mass used to form stars has been used to fuel the supermassive 
black hole. 

Another alternative would be for the mass loss by the stars 
to be the main source of fuel for the AGN. The mass loss rate for 
this population is 1 Moyr" 1 , which is larger than the mass needed 
to maintain the AGN activity, providing that the process of trans- 
ferring ejected mass from stars to the black hole is reasonably effi- 
cient. In either case, it is, nevertheless, still uncertain how this mass 
can move inwards to < 1 pc scales to fuel the black hole. 

As a cross check, we investigated the influence of consider- 
ing a larger characteristic timescale in our STARS models. We re- 
peated the above calculations for a timescale of 100 Myr. The con- 
version between the supernova rate measured and the age from the 
STARS model will be different, and will provide an upper limit for 
the age (2.5 x 10 8 yr). The mass-to-light ratio would be M/Lh = 
1.3Mq /Lq, and a mass loss rate of Mi oss = 2.1Mg/yr. The con- 
clusions remain unchanged, the mass loss rate of the stars would in 
principle be enough to fuel the black hole, or, if the gas that created 
the stars also fuelled the black hole, the fraction of gas for the AGN 
compared with gas used to form stars is again around 1 per cent. 
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Figure 12. Comparison between the symmetrised observed (left) and the JAM predicted V rms (center) for the best-fit set of parameters. The absolute values 
of the residuals are shown in the right panel: |V I ^ A i M — V r ms|- The colour-bar is in km/s. 



In summary, it is possible that the counter-rotating core is as- 
sociated with new inflow of gas which fuelled the formation of a 
new stellar population with distinct stellar kinematics. The [Fe II] 
emission traces the supernova and is observed in the inner regions 
where the counter-rotating core is. We used the measured [Fe II] 
flux and the STARS code to learn more about the episode of star 
formation and its relation with the black hole fuelling. The star- 
burst has an age of approximately 6.5 x 10 7 yr, formed 2 x 10 8 M Q 
stars and presents a mass loss rate of 1 Moyr -1 . If the AGN activ- 
ity is fuelled by ejections from the stars, we would need a mass loss 
comparable to the mass accretion rate of the black hole, which is 
0.022 M Q yr _1 . The mass loss rate obtained using STARS would 
be enough to fuel the black hole activity, providing that the pro- 
cess is efficient. On the other hand, if the new inflow of gas was 
responsible for the star formation and for providing the fuel to the 
black hole, we would need an amount of gas capable of powering 
the AGN for 6.5 x 10 7 yr. Considering that the black hole was al- 
ways active and at a constant mass accretion rate, the total initial 
gas mass needed for black hole fuelling would be 1.4 x 10 6 Mq 
which is roughly 1 per cent of the total mass used to form stars. 
These calculations provide approximate values only, but are useful 
aids to understanding the general mass budget in the central region 
of the galaxy. 



3.4 Black hole mass 

The spatial zone where the black hole potential dominates over the 
galaxy potential, the black hole sphere of influence, is given by the 



d2005h of 93.5 ± 8.5 km s . The values agree marginally within 
the errors. The difference could also be related with the distinct 
stellar population that we probe using infrared observations com- 
pared with the one observed in the optical. The method of inferring 
the black hole mass based on the velocity dispersion is known to 
be subject to large uncertainties. In the case of MCG-6-30-15, the 
decoupled kinematic components in the centre of the galaxy could 
also affect the velocity dispersion measurement. 

As a first approach, we use the velocity and velocity dispersion 
at the closest resolved region around the black hole to determine 
an upper limit for the black hole mass. The enclosed mass within 
radius R is, from virial arguments, M enc — (v 2 + 3 x a 2 )R/G. 
We determine the velocity dispersion from the integrated spectra 
in the inner R < 0".2 region around the black hole, excluding 
the R < 0" .1 due to the high non-stellar continuum. The value 
obtained is a = 89 ± 8 km/s, and the typical velocity is 10 km/s 
which gives an enclosed mass of M enc — 1.7 x 10 s Mq. From the 
calculations in Section 13.31 we can obtain a lower limit estimate 
(because it is only the mass of the young stellar population) for the 
stellar mass using the observed luminosity in the /f-band within a 
radius of R < 0".2 and the mass-to-light ratio determined from 
STARS: M st ei = 2.4 x 1O 7 M . An upper limit for the black hole 
mass can be found by subtracting the stellar mass from the enclosed 
dynamical mass: Mbh < 1-5 x 1O 8 M . 

Alternatively, we try to constrain the mass of the black hole 
and the dynamics of the galaxy with the Jeans Anisotropic Model 
(JAM) method of ICappellaril 120081 . This model generalises the 



asymmetric Jeans equations to the case of anisotropy by including 
radius from the black hole: R = GMbh/o" 2 - For MC G-6-30-15, an anisotropy parameter f3 z = 1 — erf /a,?. The model we use as- 



using the upper limit values from McH ardv et al.|j2005l) , the radius 
of influence is R ~ 3 pc or 0".02. This scale is not resolved in 
our data, since our spaxel size is 0".05 and the PSF has a FWHM 
~ 0".l. The galaxy potential will contribute significantly to the 
measured dynamics, which makes it harder to determine the black 
hole mass directly from the central stellar kinematics, without do- 
ing a multi-component dynamical model. Previously, the black hole 
mass for this galaxy was determined from the black hole relation 
with the hos t galaxy (Mbh - op. The effective radius in this galaxy 
is R e ~ 9" (Boi sson et alj2002h . lMcHardv et alH2005t) use the re- 
lation from Ferrarese (2002J) defined at a radius of R e /8 ~ 1".12 
and velocity dispersion from long slit spectroscopy to obtain a 
black hole mass of (3 - 6) x 10 6 Mq . We determine the veloc- 
ity disp ersion from a p s eudo slit with the same size as the one 
used by McHardv et "all fe005h . The integrated va lue is a = 109 
km/s, which is slightly higher than that obtained by McHard v et al.l 



sumes axisymmetric geometry and takes as input the galaxy surface 
brightness to generate a prediction for the second velocity moments 
(K- ms = \/V 2 + a 2 ), for a combination of physical parameters: 
inclination, black hole mass, mass-to-light ratio and /3. The output 
model is then compared, and adjusted if necessary, to the observed 
kinematic map of the galaxy. The galaxy surface brightness is given 
as a combination of gaussians, parametrised by a Multi-Gaussian 
Expansion (MGE) fitting software developed by Cappellaril l2002l . 
We start by doing a multi-Gaussian expansion of our /f-band map, 
constrained to a region of r < 0".8 to exclude the outer counter- 
rotating component. The surface brightness can be described by 
three gaussians, we plot the resulting contours in Fig.QT] The gaus- 
sian parameters are converted to the units required by JAM using 
the guidelines from the JAM code release by Michelle Cappellari. 
We are conservative in the number of bins we fit, and exclude the 
central pixels that show the largest errors in the velocity disper- 
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Figure 13. Compaiison between the symmetrised observed (left) and the JAM predicted Vims (centre) corresponding to a black hole mass upper limit of 
1.5 X 10 8 Mq and the remaining parameters similar to Fig. 1121 The absolute value of the residuals are shown in the right panel: |V^.^ M — V r ms|- The 
colour-bar is in km/s. 



sion determination. The errors in velocity and velocity dispersion at 
each bin were determined using the same method as in Section |3T71 
We run 100 Monte Carlo simulations for each bin, setting the error 
statistics to be the same as the data at that position. The result is 
a two-dimensional map of velocity and velocity dispersion errors 
that are taken as an optional input map for the Jeans modelling. 

The model output and parameters are affected by the area in 
which the modelling is done. A change in the black hole mass af- 
fects mainly the central bins and does not change much the Vrms 
map in the outer regions. The lower limit for the black hole mass is 
difficult to determine, since the Vrms is consistent with a null black 
hole mass (the x' 2 of the model is not very sensitive to variations 
in the mass below ~ 10 7 M o ). The black hole mass value and the 
mass limit will depend on the area where we calculate the \ 2 - F° r 
a larger area, and hence more degrees of freedom, the reduced x 2 
will increase more slowly. We use the region r < 0" A to limit the 
analysis to the area of higher black hole influence and at the same 
time to have enough bins for the calculation. The least-squares fit- 
ting routine MPFIT is used to search the parameter space for A/bh , 
inclination, j3 and M/L and determine the best fit parameters. We 
then fix all parameters, except the black hole mass, to their best fit 
values and determine an upper limit for the black hole mass based 
on the x 2 — Xmin distribution. The absolute values of the parameters 
in the fit are hard to constrain, due to possible degeneracies and the 
effect of the radius of the area analysed. The degeneracy between 
the inclination and fi z parameters for example, could be remove d 
with observationally motivated constraints (e.g. ICappellaril i2008). 
Unfortunately we do not have external information on the inclina- 
tion of the galaxy. Following the arguments presented above, we 
focus on using the model to determine an upper limit for the black 
hole mass. 

The lcr confidence limit gives an upper limit for the black hole 
mass of Mbh < 6 x 1O 7 M0. The best-fit black hole mass is 4 x 
1O 6 M0 and a comparison between the data and model is shown in 
Fig. [12] However, this best-fit value should be taken with caution 
due to the arguments discussed above; the upper limit provides a 
stronger constraint on the black hole mass. 

An increase in the black hole mass causes the Vrms of the pre- 
vious analysis to increase in the inner regions of the galaxy. In 
Fig. UH we plot an example of this effect, by fixing the parame- 
ters to be the same as in Fig.[T2]but increasing the black hole mass 
to the upper limit determined from the simple calculation with the 
integrated velocity dispersion: 1.5 x 1O 8 M0. It is clear that the in- 
crease in Vrms expected from models with high black hole mass is 
not observed in the data, which confirms qualitatively the value de- 



termined as an upper limit. This upper limit was not derived based 
on the M-cr correlation, and hence provides an independent mea- 
surement. Nevertheless, if we take the upper limit for the black hole 
mass (6 x 1O 7 M0) and the integrated velocity dispersion (a — 109 
km/s) we measure, the location of MCG-6-30-15 i n the Mbh — o 
plot i s in agreement with the Mbh — o relation (e.g. lGiiltekin et all 
2009). We conclude that the mass of the black hole in MCG-6- 
30-15 is lower than 6 x 1O 7 M0 which is in agreement with the 
previous estimates. 

For a bolometric luminosity of 1.3 x 10 44 erg s _1 
dVasudevan et alj|2009h . we obtain a lower limit for the Eddington 
ratio of A > 0.02. 



4 CONCLUSIONS 

In this work we studied for the first time the inner ~ 470 pc of the 
galaxy MCG-6-30-15 using integral field spectroscopy. We were 
able to remove the AGN broad hydrogen Brackett emission lines 
which were dominating the spectra in the central region of the field- 
of-view, and measure the properties of the stellar absorption lines. 
The stellar kinematics of this galaxy can be characterised by a low 
rotational velocity (~ 40 km/s) compared with the velocity disper- 
sion (80 — 140 km/s). The velocity dispersion is higher at larger 
radii and close to the major axis and decreases when approaching 
the position of the black hole. There is a change in the direction 
of the stellar rotation when comparing the central core (r < 0".8) 
with the outer regions. We argue that we are in the presence of a 
galaxy with a counter-rotating core, due to the rotation curve and 
the observed drop in velocity dispersion at small radii. The gas dy- 
namics traced by the [Fe II] emission line show an asymmetric dis- 
tribution in the inner r < 0".8 arcsec, with a higher intensity and 
larger velocity dispersion on the West side of the nucleus. It also 
shows a velocity gradient with blueshifted velocities of —30 km/s 
on the West side and +120 km/s on the East side of the nucleus. 

The quality of the data does not allow us to exclude the pos- 
sibility of an outflow, but due to the smooth velocity curve and its 
spatial distribution we argue that the [Fe II] has been excited by 
supernova shocks. In this scenario we use a model to reproduce the 
supernova rate inferred and determine the star formation history. If 
the counter-rotating core is a result of a recent inflow of gas that 
formed the stars and led to the supernova explosions, we can deter- 
mine how much gas was used to form stars and how much gas was 
used to fuel the black hole. The percentage of gas used to fuel the 
black hole is at most 1 per cent of the gas used to form stars. If on 



© 0000 RAS, MNRAS 000, 000-000 



Black Hole and Stellar Population of MCG-6-30-15 1 3 



the other hand the outflows from stars fuel the black hole, we con- 
clude that with this supernova rate, the winds from stars would be 
enough to fuel the AGN during the age of the starburst ~ 6.5 x 10 7 
yr. These arguments are of course dependent on how efficient the 
process of transferring gas from larger scales to the black hole is, 
but can give us a general overview of the mass budget in the vicinity 
of the nucleus. 

Using the measured kinematics at r < 0".2, we are able 
to determine an independent upper limit for the black hole mass 
of 1.5 x 10 8 Mq, which is consistent with other estimates. We 
also reproduce our observations using a dynamical model, and 
determine an upper limit based on the model predictions for the 
Kms = VV 2 + a 2 in the central r < 0".4 of the galaxy, of 
M BH < 6 x 10 7 M Q . 

The study of MCG-6-30-15 allowed us to determine the dy- 
namical properties of the inner regions of this galaxy and infer its 
stellar history. There are a growing number of galaxies that have 
been observed in detail using integral field spectroscopy. MCG-6- 
30-15 has a larger mass accretion rate (high AGN activity), which 
is not common among the usually selected targets due to the dif- 
ficulty in removing the AGN contamination. We have shown here 
that it is possible to obtain the stellar kinematics with the presence 
of Brackett broad emission lines, allowing in the future to increase 
the parameter range in AGN activity of the galaxies studied. The 
kinematically distinct core in the centre of the galaxy may be asso- 
ciated with bar-driven gas inflow, which could, on a smaller scale, 
be related with the fuelling necessary for AGN activity. We relate 
the available gas to form stars with the observed AGN activity and 
obtain general constraints on the mass budget in the centre of the 
galaxy. Combined studies of stellar properties and AGN activity in 
the centre of galaxies, will in the future help clarifying the relation 
between star formation and black hole fuelling at small scales. 
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